ON THE WELL-POSEDNESS OF A MATHEMATICAL MODEL 
DESCRIBING WATER-MUD INTERACTION 
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Abstract. In this paper we consider a mathematical model describing the two-phase 
interaction between water and mud in a water canal when the width of the canal is small 
compared to its depth. The mud is treated as a non-Netwonian fluid and the interface 
between the mud and fluid is allowed to move under the influence of gravity and surface 
tension. 

We reduce the mathematical formulation, for small boundary and initial data, to a 
fully nonlocal and nonlinear problem and prove its local well-posedness by using abstract 
parabolic theory. 



1. Introduction 

We consider herein a mathematical model describing the evolution of a two-phase system 
consisting of mud and water. This problem is of immense physical importance as it is 
related to sediment transportation in coastal and estuarine regions. Therefore, a large 
number of experiments and theoretical investigations have been and are still dedicated to 
better understand two-phase mud-water systems. 

We shall treat the mud as being a non-Newtonian fluid, which reflects the rheological 
behaviour of fluid muds, cf. [19], as the viscosity changes with the rate of strain. The 
viscosity function for the mud, though general, is assumed to satisfying certain, not to 
restrictive, conditions. These conditions appear in the literature in the context of the non- 
Newtonian Navier-Stokes problem [2] as well for the Stokesian Hele-Shaw flow [6] and it 
was proven that these are sufficient conditions to ensure the local well-posedness of these 
problems. A diffuse interface model for the flow of two viscous incompressible Newtonian 
fluids in a bounded domain have been analyzed at the level of weak solutions in [1], when 
considering fluids with different densities. Assuming a sharp interface between the fluids, 
which are modeled by the full Navier-Stokes equations and have constant densities, the 
new formulation has been studied in [14, 15] by using maximal regularity for the linearized 
problems. 

Numerically, this physical setting has also received plenty of attention leading to a large 
number of models and experiments investigating the entrainment and settling of the fluid 
mud, cf. [9, 11. 20]. Other experimental studies are concerned with the damping effect of 
fluid muds on the waves traveling at the water surface, see e.g. [4, 21]. 

Fluid mud is a highly concentrated suspension of fine grained sediments and is treated 
herein as being a homogeneous fluid with general shear dependent viscosity. Furthermore, 
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the width of the canal is assumed small compared to its length and hight, so that, width 
averaging the Navier-Stokes equations we can approximate the flow in the water and mud 
region by two-dimensional linear and nonlinear Darcy's laws, respectively, cf. [10]. In 
this way we incorporate the viscous behaviour of the fluid mud into our model but loose 
the generality one has when working with the full Newtonian and non-Newtonian Navier- 
Stokes equations. Assuming also that there is a sharp interface between mud and water, 
we are lead to a coupled problem describing the evolution of the water-mud system, and, 
if the mud is replaced by a Newtonian fluid, we rediscover the Muskat problem [3, 7, 22] 
which is a current research topic. We note also that there are recent investigations [13] of 
the least squares formulation for a two-phase coupled problem with Stokes flow in the one 
subdomain and a linear Darcy flow model for the second fluid. 

Having a general viscosity function makes the analysis more involved compare to the 
Newtonian case [3, 7, 22]. Using elliptic theory for quasilinear equations though, we are able 
to formulate the problem with the help of a nonlinear and nonlocal operator depending 
on the function parameterizing the sharp interface between the mud and water and its 
time derivative. Taking advantage of the properties of the flow, that is the volume of 
fluids are preserved, the implicit function theorem allows us to reformulate the problem 
as an abstract evolution equation for small boundary and initial data. Under a smallness 
condition for the time derivative of this function, which is necessary only if the mud is a 
shear thinning fluid, we prove then by using parabolic theory [16] existence and uniqueness 
of local solutions for our problem. 

The outline of the paper is as follows. In Section 2 we present the mathematical model 
and state the main result of this paper, Theorem 2.1. In Section 3 we transform the coupled 
system of equations obtain in Section 2 and express the problem as a single operator 
equation. The proof of the main result is found in Section 4. 



2. The mathematical model and the main result 

We present now the mathematical model describing the two-phase interaction between 
water and mud in a water canal with small width. The subscript m is used when referring 
to the mud and which occupies the domain Q m (t), respectively w for the water located in 
£l w (t). These two fluids are separated by the interface T(t) which evolves in time and is to 
be determined as a part of the problem. The physical relevant problem is essentially three- 
dimensional, but it may be approximated by a two-dimensional mathematical model. We 
refer to [10] for a deduction of a generalized Darcy's law for non-Newtonian fluids, which 
serves as a approximation of the non-Newtonian Navier-Stokes equations. Thus, we restrict 
our considerations to the situation when Q m (t),£l w (t) C M 2 . 

The mud-water problem is a potential flow, in the sense that the velocity fields v w 
satisfy Darcy laws. The motion of the mud is governed by the nonlinear Darcy's law 



Vu r , 



/«m(|Vli m | 2 ) 

and in the water domain we set 



in n m (t) (2.1) 



-— Vu w in Q w (t), (2.2) 
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fj, w being the constant viscosity of the water. The effective viscosity \i m is defined (cf. [10]) 

by 

^(P):=y_^ctaj forr>0. (2.3) 

With /i € C°°([0, oo), (0, oo)) denoting the shear-rate dependent viscosity of the mud, we 
let jl : = c/i o [r i— > r / u 2 (r)]~ 1 , whereby c is a constant depending on the gap width of the 
canal. The invertibility of the mapping [[0, oo) 3rn> rp?{r) £ [0, oo] is guaranteed if we 
restrict our considerations to fluids having viscosity functions which satisfy: 

m< n< M and m < fi(r) + 2r//(r) < M (2.4) 

for some positive constants m, M. It follows from (2.3) that fx m inherits similar properties 
as n, that is jjL m € C°°([0, oo), (0, oo)) and \x m satisfies 

m < Mm < M and m < /i m (r) — 2r/i m (r) < M (2-5) 

possibly with different constants. We refer to [6] for precise calculations. Moreover, the 
scalar functions u w and defined by 

u m ■= Pm + gp m y 'm. Q m (t) and u w :=p w + gp w y inU w (t), (2.6) 

whereby p m ,p w are the dynamic pressures, p m , p w denote the constant densities of the 
fluids, g is the gravity constant, and y is the height coordinate. Relations (2.4) are satisfied, 
besides by all Newtonian fluids, by the Bingham model with dynamic viscosity 

Kr) = Moo + ^ (2.7) 

if and only if (3tq < ^jj,^. For this reason, suspensions of hectorite can be considered in 
our model (cf. [19]). 

Additionally to (2.6), we assume that the fluids are incompressible, that means 

V • v m = in Q m and V ■ v w = in Q w . (2-8) 

The surface tension plays a major role in interfacial phenomena as a restoring force. It 
compensates the pressure difference across the interface T(t) according to Laplace- Young's 
condition, so that we get the following equation at the interface between the fluids: 

u w ~ u m = 7«r(t) + g{Pw - Pm)y on F(t), (2.9) 

with 7 being the surface tension coefficient at the interface, and «r(i) i s ^ ne curvature of 
r(i). Equalizing the normal components of the two velocity fields on the interface T(t), we 
also obtain two kinematic boundary conditions at this boundary: 

V(t,-) = (v m (t,-),u(t,-)) = (v w (t,-),u(t,-)) on r(i), (2.10) 

which means that the interface moves along with the fluids (the interface consists of the 
same particles at all times). We denoted here by u(t) the unit normal at T(t) which points 
into Q w (t), and V(t) is the normal velocity of the interface. To complete the model, we 
presuppose that the bottom r_i = [y = — 1] of the canal is impermeable, that is 

u m ,2 =0 on r_ 1; (2-11) 
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where u m ^ and u mt ij, 1 < i,j < 2, stand the first and second order derivatives of u m 
(similarly for u w ). Also, we prescribe Dirichlet boundary conditions at the interface Ti := 
[y = 1], which is presupposed to be located above T(t), that is: 

u w = h(t,x) on F±. (2-12) 

Lastly, we assume the interface at time t = to be known 

r(o) = r . (2.13) 

System (2.1)-(2.13) is a two-phase moving boundary problem. The main interest is to 
determine the motion of the interface V{t) separating the fluids. If we know T(t), then we 
can find the potentials u m and u w by solving elliptic mixed boundary value problems. 



2.1. Parameterizing the boundary. In the following we restrict our considerations to 
horizontally 2ir— periodic flows. The unknown interface T(t), which is to be determined, is 
assumed to be parametrized by a function from the set 

W:={/eC(S):||/|| c(S) <l/2}, 

where S = R/2tt7j is the unit circle. Let a £ (0, 1) be fixed for the remainder of this 
paper and define X := /iq +q (§) when we incorporate surface tension effects into our prob- 
lem, respectively X := /jq +q (§) when 7 = 0. Moreover, we let Y := /iq + "(§). The small 
Holder space /i m+/3 (§), m G N and ft G (0, 1), is defined as being the closure of C°°(§) in 
C m+ "(§). Particularly, small Holder spaces are densely and, by Arzela-Ascoli's theorem, 
also compactly embedded in those with lower exponent. A further important feature in 
our context is that they are stable under the continuous interpolation functor (•|-)2 00 , cf. 

[16]. The subspaces /i™ +/3 (§) consist only of the function in /i m+ ^(§>) with integral mean 
zero. Our choice for X and Y is due to the following observation: if T(t) is at rest the 
flat line y = 0, because there is no flux over cf. (2.11), the volume of fluid mud of a 
solution T(t) should match the one at rest. Lastly, we define V := U (~l X and observe that 
is an open neighborhood of the zero function in X. 

By our modeling considerations, the two fluids completely fill the domain fi := S x 
( — 1,1). If / : [0,T] — s> V, with T > 0, is a function describing the evolution of the 
interface between the fluids, then at any time t € [0,T] we have that £l m (t) = O m (/(t)), 
tt w (t) = tt w (f(t)), and r(t) = r(/(t)), where, given h € V, T(h) := {(x, h(x) : x G §} and 

Mm(h) := {(x,y) GO: -1 < y < h(x)}, n w (h) := {(x,y) G O : h(x) < y < 1}. 

In order to reformulate the kinematic conditions (2.10) we follow the evolution of a particle 
(x(t),y(t)) on the interface r(/(£)). By the definition of T(f(t)), we know that y(t) = 
f(t,x(t)), and differentiating this relation with respect to time yields Off = ( — f',l)v, 
where v(t) is the velocity of T(t). This shows that the normal velocity V(t) of the interface 
r(t) is given by V(t) = d t f{t)/yjl + f' 2 (t). 
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Therewith, the evolution of the water-mud system (2.1)-(2.13) is described by the fol- 
lowing system of equations: 

Au w = in 0, w (f), 

v( , V " m |2 = in n m (f), 

u w = h on Ti, 

u m ,2 = on (2.14) 

u w -u m = 7«r(/ ) + g(pw - p m )f on r(/), 

dtf = V 1 + f' 2 dv u w = ^77^ 7^d u u m on 

/(0) = /o, 

where /q G V is chosen such To = T(/o). We shall assume that the Dirichlet boundary data 
h is continuous in time, more precisely 

h G C([0,oo),/i 2+a (S)). (2.15) 

For the sake of simplicity, we identify functions defined on T± and T_i with functions on 
the unit circle S. 

We shall call a triple (f,u w ,u m ) classical Holder solution of (2.14) on [0,T] if 
/GC([0,T],V)nC 1 ([0,T],Y), 

u w (t,-) € h 2+a (n w (f(t))), u m (t,-) G /i 2+a (a m (/(t))), i G [0,T], 

and if (f,u w ,u m ) satisfies (2.14) pointwise. Given functions /i,/2 G C(§) with /i < /2> 
we let Qf 1 j 2 '■= {{x,y) '■ fi(x) < y < f2(x)}. The space h 2+a {Q,f 1 j 2 ) is defined then as 
being the closure of the smooth functions C°°(^/i,/ 2 ) m the Banach space C 2+a {Q,f 1 j 2 ). 
The main result of this paper is the following theorem: 

Theorem 2.1 (Local well-posedness). Let 7 G [0,oo), li£K and assume that 

7 > or p m > p w . (2.16) 

Then, there exist Of an open neighborhood of zero in X, Z an open neighborhood of zero in 
Y, and O-^ an open neighborhood of h in /i 2+a (§) such that if fo G Of, h satisfies (2.15), 
and h(0) G O-^ then there exists T > and a unique maximal solution f : [0, T) — > Of 
satisfying dtf{t) G Z and h(t) G Oj^ for all t G [0, T). Moreover, if p! m > then Z = Y. 

Remark 2.2. In contrast to the constant viscosity case studied in [7] where, at any time 
to where the solution exists, the time derivative dtf (to) is uniquely determined by the 
boundary conditions and f(to), the situation considered herein is more complex and it is 
not clear whether this feature is preserved when the viscosity fi m is decreasing. Particularly, 
when considering the fluid mud model (2.7) from [19], dtf is uniquely determined by / 
under the smallness condition dtf G Z, cf. Proposition 4.2. 

3. The transformed problem 

The difficulty when dealing with problem (2.14) is mainly due to the fact that the equa- 
tions of (2.14) are defined on domains which evolve in time. To overcome this problem, 
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we shall transform the system by writing its equations on fixed domains. This will trans- 
form though the differential operators into operators with coefficients depending on the 
unknown /. Then, we will introduce solution operators to linear and quasilinear elliptic 
boundary value problems involving these transformed differential operators and study their 
regularity properties. At the end we obtain at an operator equation for the function / and 
its time derivative. 

To simplify our notation, W6 introduce first th.6 cjUctsilinecLr opera/tor Q \ C^(Q m ) — y 
C(Q m ) by setting 

Qu m : = V ( ,Y^ m l9 . ) for u m € C 2 (S1 



/"m(|Vtl m | 2 ) 



ni ) 



For u m e c 2 (n m ), one easily computes that Qu — Or/(Vu m )it mj jj, with coefficients Gjj, 1 < 
i, 3 < 2, given by 

Given z £ M 2 , the eigenvalues of the matrix (aij(z))i<ij<2 are 

a ( Z ) = 1 a ( Z ) = 1 _ 2 N 2 ^(N 2 ) 
M{z) » m (\z\ 2 r [ ) MM 2 ) /4(M 2 ) ' 

and it follows readily from the fact that \i m satisfies the relations (2.5) that the eigenvalues 
are positive and uniformly bounded in z £ M 2 , i.e. Q is a uniformly elliptic quasilinear 
operator. 

From now on, we set £l m : = O m (0), £l w := Q w (0), and will identify the common bound- 
ary To := S x {0} of fl m and fl w with the unit circle S. Given / £ V, we define the mapping 
<f>t : ft — > fl by the relation 

4>f{x, y) ■= (x, y + (l- y 2 )f{x)) for (x, y) £ O. 
Since functions / £ V satisfy ||/||c(s) < 1/2, one can easily verify that the function 



is the inverse of (j)f , that is (j>f is a diffeomorphism having the same regularity as /. More- 
over, we have that (f)f(J} m ) = £l m (f), 4>f(Cl w ) = fl w (f), an d <^/(Fo) = To keep the 
notation simple, we shall also denote by 4>f t ne diffeomorphisms obtained by restricting (ftf 
to Q m and Q w . 

Using these diffeomorphisms, we re-express now the equations of (2.14) on the fixed 
reference domains Q m , Q w , and their boundaries. To this end, we introduce for each / £ V 
the operator A w (f) : h 2+a (fl w ) — > h a (fl w ) by setting 

A w (f)v w := (A(v w o (pj 1 )) o <j)f for v w £ h 2+a (fl w ). 
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A simple computation yields that 

AvU)v w =-5-5- + 2- + ~~ 7i OT77T2 - + 



dx 2 1 - 2y/ dxdy \ (I - 2yf) 2 (1 - 2y/) 2 J dy 
7V-1) + 4/' 2 y(y 2 - 1) , 2//' 2 (y 2 - l) 2 + 2/^ 



1 - 2yf (1 - 2y/)2 (1 - 2yff J dy ' 

for v w S /i 2+a (O u ,). Consequently, one can see that A w (f) is uniformly elliptic and that 
it depends analytically on /, that is A w € C W (V, C(h 2+a (Q w ), h a (Q w ))). Furthermore, 
corresponding to the operator Q, we let A m : V x h 2+a (Q m ) — > h a ({l m ) be the operator 
with 

A m (f,V m ) := \Q(v m ° ° <Af for Vm G ^ 2+Q (^m)- 

Lemma 3.1. The operator A m is smooth Am G C°°(V x /i 2+a (r2 m ), h a (Q m )) and has a 
quasilinear structure A m (f, v m ) = bij(y, f,VfV m )v m ^j + b(y, f, V/f m )t» m>2 whereby 

( f'(y 2 -l) 1 
V/U m := I f m ,l + 1 _ 2y / 1 _ 2yf Vm ' 2 

and the coefficients of the operator are given by 

b u (y,f, Vfv m ) =au(Vfv m ), 

f'( y 2 - 1) 1 
bi2{y,f,VfV m ) = 1 _ 2y f Q n( V /^) + I _ 2y 7 ai2 (y ' f v m), 

&22(y, /, V/« ro ) = ( 1 _ 2y /)2 °H( V /^) + 2 (1 _ 2y/) 2 Q 12( V /^m) + 1 _ 2yf a22 ( V f l 



b(y, f, VfV r , 



f"(y 2 -l) , 4/' 2 y(y 2 -l) , 2/f 2 (y 2 -l) 2 ^ 



+ 7Z ' + 



an(V/w m ) 



1-22// (l-2y/) 2 (l-2yfY 
4f'y-4ff'(l + y 2 ) 2f 
+ (1 - 2yfr ai2(V ^ m) + (l- 2y /)3 a22(V ^ ) ' 



forfev, v m eh 2+a (n m ). 

Proof. The quasilinear structure of A m is found by using the definition of Q and A m . Then, 
taking into account that /i m is smooth, we deduce that A m is a smooth function of both 
variables. □ 

Corresponding to the fifth equation of system (2.14), we define the boundary operators 
B w (f ) : h 2+a (Q w ) -> h 1+a {S) and B m : V x /i 2+Q (O m ) -> /i 1+a (S) by setting 

s«(/K = — tr ((vk o !(-/', 1)) o <t>/j , Vw g h 2+a (n w ), 



-1- 



S m (/,« ro )=tr (—77^7 T^TT^-^ 1 ) W/ h (/,geVxft 2+a ((! m ), 

W/imdV^mO^ )| 2 ) / J 
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respectively, where tr is the trace operator with respect to § = Tq. More precisely, we have 
that 

B w (f)v w =fj, w 1 [(l + f' 2 )tTv w , 2 -f'tvv w , 1 ] for v w € h 2+a {VL w ), (3.1) 

so that B w (f) depends analytically on /, i.e. B w € C U (V, C(h 2+a (n w ), h 1+a {S))). Con- 
cerning the boundary operator B m , we get 

B m (f,v m ) = {1 + f " ) ° m l*~ f 'n Vm ' 1 for (f,v m ) S V x h 2 + a (Q m ), (3.2) 

Hm (tr \\/ fV m \ ) 

whereby trV/i> = (v m ,i — /V,2,%,2). From (2.5) together with the smoothness of fj, m , 
we deduce that B m is a smooth operator B m G C°°(V x h 2+a (£l m ) , h 1+a (§)) . Letting 
v w = u w o and u m = u m o<f>f, we see that the tuple (f,v w ,v m ) satisfies the following 
system of equations 



dtf 





= 


in 




■Am\f ■> 


= 


in 




V w 


= h(t,x) 


on 


Ti, 




= 


on 


r_i 


Vw V m 


= 7K(f)+9(Pw-Pm)f 


on 


r , 


-B w (f)v w 


= -B m (f,v m ) 


on 


r , 


/(0) 


= fo, 







(3.3) 



pointwise. When 7 > 0, we set k : V -> /iQ +a (§) to be defined by «(/) = /"/(l + / /2 ) 3/2 
for / € V. Let us observe that the curvature operator k is also a real- analytic function 
K € C^(V,/io +Q (S)). We enhance that the equations of (3.3) are all defined on sets which 
are fixed in time. 

Defining the notion of classical solutions for (3.3) in a similarly way to that for (2.14), 
we have the following equivalence between problems (2.14) and (3.3). 

Lemma 3.2. If (f,u w ,u m ) is a solution of (2.14), then the tuple (f,u w o(f)f,u m o(f)f) is a 
solution of (3.3). Conversely, if(f,v w ,v m ) is a solution of (3.3), then (/, u^o^T 1 , v m o(f)J ) 
is a solution of (2.14). 

Proof. Assume that (f,u w ,u m ) is a solution of (2.14). As we already noticed, the tuple 
(f, u w 4 > fi u m <fif) solves (3.3) pointwise. We show now that v w := u w (f>f has the 
required regularity v w € h 2+a (£l w )- To this end, let (/„) C C°°(§) be a sequence such that 
/ n — >■ / in C 2+a (§). We can choose (/„) such that in fact f n > f for all n £ N (we do not 
require that f n € V). If (u w ) £ C°°(£l w (f)) is a sequence such that u w — » u w , then due 
to £l w (f n ) C £l w (f)i the function u w <fif n is well-defined and u w 4>f n £ C°°(r2 u ,). The 
function 7;^ is an accumulation point of the set {u w 4>f„ '■ P,n £ N} : 

IK <t>fn ~ «wlloa+«(n^ ^ IK ° ~ <> ^/Hc2+-(n^) + II « ~ u ™) ° 0/llc 2 +"(?^) 

< C (ll n S>llc 4 (n^)ll/« ~ /llc 2 +-(§) + +\\ u w ~ u w\\c^(Tu)) ' 

so that indeed v w € h 2+a (£l w ). Analogously, v m £ h 2+a (Q m ) and we proved the first claim. 
The reciprocal implication follows by similar reasons. □ 
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3.1. The abstract formulation of the problem. We enhance that even if the transfor- 
mation performed in the previous paragraph has the hindrance of introducing additional 
nonlinear coefficients, it allows us to re-write the problem as a single nonlinear equation. 
To do this, we define in the following solution operators to elliptic boundary value problems 
which are closely related to system (3.3). They are used later on to formulate the original 
problem (3.3) as a single operator equation. 

Lemma 3.3. Given f G V and (F, h) G h 1+a (S) x h 2+a (8), the problem 

A w {f)v w = in Q w , 

B w {f)v w = F on r , (3.4) 

v w = h on T\. 

possess a unique solution T(f)[F, h] G h 2+a (£l w ). Moreover, the operator T is real-analytic 
T G C^(V,£(/i 1+Q (S) x h 2+a {§),h 2+a (n w ))). 

Proof. Given / G V, the problem (3.4) possesses a unique classical solution u w o(/)f, where 
u w is the unique solution of the linear elliptic boundary value problem 

Atta, = in Vl w (f), 

^s/l + f^d^ = F on T(/), (3.5) 



u w = h on T[. 

Denoting by tri the trace operator with respect to T±, the operator T(f) is the inverse of 
the linear operator 

(A,(/),-M/),tri) : C 2+a (T^) C a {TK,) x C7 1+a (S) x C 2+Q (S). 

Since A w and B w depend analytically of /, and the function associating to a bijective oper- 
ator its inverse is also real-analytic, we get T G C"{V, £(C 1+a (S) x C 2+a {§) , C 2+a {TQ)). 
Due to elliptic regularity, cf. Theorem in [8], the solution v w of (3.4) belongs to C°°(£l w ) 
if /, F, h are smooth functions. By continuity, we conclude that T(f)[F, h] € h 2+a (VL w ) for 
all / G V and (F, h) G h 1+a (§) x h 2+a (§), which completes the proof. □ 

We enhance that second equation of (3.4) corresponds to one of the kinematic boundary 
condition in (2.14). Whence, if we know / and its derivative dtf, the function v w is given 
by v w = T(f)[—dtf,h]. We can use now v w as a boundary data in the fifth equation of 
(3.3) and determine v m as a function of / and dtf. 

Lemma 3.4. Given f G V and p G h 2+a (S), the problem 

AmifjVm) — in 

v m = p on T , (3.6) 

v m ,2 = on r_i 

possess a unique solution 7t(f,p) G h 2+a (Q m ). Moreover, the operator 1Z is smooth 1Z G 
C°°(V x h 2+a (§),h 2+a (n w )). 

Proof. We pick / G V and consider the quasilinear Dirichlet problem 

Qu m = in 0_ 2 _/j, . . 

u m = P on 9fi_ 2 -/,/. 
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Theorem 8.3 in [12] together with condition (2.5) guarantees the solvability of (3.7) in the 
Banach space C 2+a (f2_2-/,/)- Furthermore, the uniqueness of the classical solution of (3.7) 
is a consequence of the weak elliptic maximum principle for quasilinear elliptic equations, cf. 
Theorem 10.1 in [8]. Denoting the unique solution of (3.7) by u m , because of the symmetry 
of the domain 0_2— / / an d of the boundary data, the weak elliptic maximum principle 
ensures ensures also that u(x,y) = u(x,—2 — y) for all (x,y) G 0_2_//. Particularly, 
u m) 2 = on r_i, meaning that v m : = ("U m |f2 m (/)) o (fij 1 G C 2+Q? (fl m ) is a solution of 
(3.7), which is unique by the same result, Theorem 10.1 in [8]. . Whence, for all / € V 
andp G h 2+a (§), the problem (3.4) possesses a unique solution TZ(f,p) G C 2+a {Vt m ). Since 
T^ifiP) £ C°°(0, m ) when both functions / and p are both smooth, if we show that 7Z is a 
smooth operator, then, by a density argument, we obtain that 7Z maps into h 2+a (Q m ). 

In order to prove this regularity property, we notice that (A m (f, v m ),tT v m —p, c^m) = 
if and only if v m := lZ(f,p), where we consider 

(f,P,v m ) -> (A m (f,v m ),trv m -p,v m , 2 ) G C a (TU) x C 2+Q (S) x C 1+a (S). 

as a map defined on V x C 2+Q (§) x C 2+a (Q m ). Since A m is smooth, if we prove that 
d Vm (A m (f, -),tr — p, d2)\ v is an isomorphism, then by the implicit function theorem 

we deduce that 1Z is a smooth map and we are done. From Lemma 3.1 we obtain that 
d Vm (Am(f,-),tr-p,d 2 )\ Vm = (T,tr,d 2 ) : C 2+Q 0U) -> C a (H^) x C 2 +«(S) x C^S), 
where 

Tw m :=bij(y, /, V fV m )w m ,ij + /, V/<; m )w m ,2 + v m ^d 3 b ij {y, f, V f v m )V f w m 
+ v m>2 d 3 b(y, f, V f v m )V f w m . 

Hereby, S3 denotes differentiation of the coefficient functions with respect to the third 
variable V fV m Since T is a linear elliptic operator and the coefficient of w m is zero, we 
conclude that (T,tr, 82) is an isomorphism, and the proof is completed. □ 

With this notation, we remark that (f,v w ,v m ) is a solution to (3.3) if and only if 
/(0) = /o, ««, = T(f)[-d t f,h], v m = K(f, trv w - jn(f) + g{p m - Pw)f), and 

dtf + B m (f,TZ(f,tTT(f)[-dtf,h] - 7 «(/) + <?(p m - *,,)/)) = (3.8) 

for all t G [0, T], for some T > 0. We shall now justify our choice of working with functions 
having integral mean zero. 

Lemma 3.5 (The volumes of both fluids are preserved). Given (/, F) G /iQ +a (§) x /iq +Q! (S), 
we have that 

F + B m (f, TZ(f, tr T(f)[-F, h] - 7 «(/) + ff (p m - Pw )f)) G ^ +Q (S). 

Proo/. Let v w := trT(f)[-F,h], v m := 1Z(f,trv w -jK(f)+g(p m -p w )f), u w : = u^o^ 1 
and u m := v m o 07 . Having assumed (2.11), we find from Stokes' formula that 

B m (f,v m )dx= [ ( ™ m \(-f',l)) (x,f(x))dx = ( d u ( Y"" 1 ) da 



V ( ^ m |2 , ) 

O m (/) V/^m(|VM m ^), 
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with da denoting the curve integral. This and our assumption F G h,Q +a (E>) leads us to 
the desired conclusion. □ 

4. Proof of the main result 

Resuming the last section, we have reduced the original problem (2.14) to the following 
problem: given /o € V find a positive time T = T(/o) > and solution curve / with 
/€C([0,T),V)nC 1 ([0,T),Y) and 

T(h(t)J(t),d t f(t)) = for allie [0,T). (4.1) 

Hereby, T : /i 2+a (S) x V x Y — > Y is the nonlocal operator given by 

T(h, f, F):=F + B m (f, K(f, tr T(f) [-F, h] - 7 k(/) + g(p m - p w )f)). (4.2) 

We note that the operator (4.2) is fully nonlinear in each of its variables. 

First, we observe that the tuple (h, 0, 0) is a solution of the equation F(h, f,F) = for 
any constant /i£l. Using the implicit function theorem we may express F as a function 
of (h, f) at least in a neighborhood of the point (h,0, 0), cf. Lemma 4.1. Doing this, we 
arrive at a fully nonlinear abstract evolution equation. Indeed, if h € K, and (/, F) = (0, 0) 
then problem (3.4) has the constant solution T(f)[— F, h] = h, and therefore we also have 
TZ(f, tr T(/) [— F, h] — 7«(/) + g(pm — Pw)f) = h. Whence, (/i, 0, 0) is indeed a solution of 
F(h,f,F)=0. 

Lemma 4.1. Given /iSR, the Frechet derivative dpj-{h, 0, 0) is an isomorphism, that is 

d F F(h,0,0) G Isom(Y). 

Proof. Using the chain rule and the linearity of the solution operator T with respect to 
(F, h), we find that 

d F T(K,o,o) = id ¥ - d Vm B m (o,h)[dpno,h)[^no)[-,o}}}- (4-3) 

It readily follows from (3.2) that we have 

d Vm B m (0,h)[w m ] = for all w m € h 2+a (fl m ). (4.4) 

Further on, using the definition of the operator 1Z and Lemma 3.1 we get that w m : = 
d p lZ(0, h)[q] is the solution of the problem 

Aw m = in £l m 

w m = q on r (4.5) 

w m ,2 = on r_i 

for all q G h 2+a (S). Pick now q G h 2+a (S) and consider its Fourier expansion q = 
Y.k&^ e%kx - Expanding also w m as w m = d p TZ(0,h)[q] = Y^k& B k{y) e%kx we b Y 
inserting both expressions into (4.5) and comparing the coefficients of e lkx for every k, 
that Bk is the solution of 

f B k (y)" -k 2 B k {y) = -1< y < 

B k (0) = b k 
B' k {~l) = 
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and therefore 

Bkiv) = + ^=k-- ky ) kez. (4.6) 

Next, we consider the Fourier expansion of F = ^feez\{o} a k e%kx G Y. Similarly as above, 
we obtain from the definition of the operator T that 

trT(0)[F,0] = Y, -^ tj ^W ikX - (4-7) 

k£Z\{0} 

Gathering the relations (4.6) and (4.7), we obtain from (4.3) and (4.4) that 

d F T(h,0,0)[F] = U + -^ta,nh 2 (k)]a k e ikx for F = ^a fc e ifca; . (4.8) 

mGZ\{0} V AW J / fcgZ 

Whence, we have shown that (9^J-"(/i, 0, 0) is a Fourier multiplier with symbol m(k) := 
1 + (/^«/Vm(0)) tanh 2 (fc), k G Z \ {0}. Since, sup fceZ \l/m(k)\ < oo, we infer from [17, 
Theorem 3.6.3] that the Fourier multiplier with symbol 1/m belongs to C(C p+a (E>)) for all 
p G N. From this, it follows at once that it belongs also to C(Y) and that it is the inverse 
of d F T(h,0,0). □ 

Using Lemma 4.1, we obtain the following reformulation of the equation J-(h, /, F) = 0. 

Proposition 4.2. Given /iGK, there exist an open neighborhood WxOx2c h 2+a (S) x 
V x Y of (h, 0, 0) and a smooth operator $:Wx0C h 2+a (S) xX->¥ with the property 
that if (h, f , F) £ W x O x Z is a solution of the equation J-(h, f,F) = then 

F = <5>{h,f). (4.9) 

Moreover, if /j>' m > 0, then we can choose Z = Y. 

Proof Since T{h,0,0) = 0, T G C°°(/i 2+a (§) x V x Y,Y), and d F T(h,0,0) G Isom(Y), 
cf. Lemma 4.1, we obtain from the implicit function theorem that there exist an open 
neighborhood WxOxZ C h 2+a (8) x V x Y of (h, 0, 0) and $ € C°°(0 x W, Z) such that 

J"(/, h,F) = and (f,h,F)eOxWxZ implies F = $(/, /i). 

In the rest of the proof, we show that for shear thickening fluids, we have Z = Y. To this 
end, let / £ V, /i £ h 2+a {§), and Fi,F 2 G Y be given such that 

T(h,f,F 1 )=T(f,h,F 2 ) = 0. 

If 

u i m = T(f)[-F i ,h)]o ( j ) j\ i = l,2 

< = ft(/, tr T(/) [-F ia fc] - 7 «(/) + 3(Pm - o i = 1, 2 
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then, by the definition of T,7Z, B w , 
system of equations 



and B m , we find that (u l m ,u l w ) solve the following 



V 





= 


iii 


««,(/) 


r>) 


= 


in 


0m(/) 


< 


= h 


on 


Tl, 


u m,2 


= 


on 


r-i, 



(4.10) 



it 



*i = \A + f' 2 d, 



Defining U„ 



m 

i 



7«r(/ ) + g{Pw ~ Pm)f 



2 d u u m 



on 

on r(/). 



u 2 and £/,, 



w 2 we obtain by subtracting the corresponding 



equations of (4.10) for i = 1 and i = 2 that 



A£/ w = in 

£f7 m = in J2 m (/), 

C/^ = on Ti, 

C/ mi2 = on r_i, 

C^-C/m = on T(/), 



(4.11) 



where 



Qu m - Qu Ti 



2 )u 2 •• 



and Cjj := J Q daij(t'S7u} m + (1 — i)Vit^) (it, 1 < i,j < 2. Let us assume that U w is not 
vanishing and its maximum is positive (the case when U w possesses a negative minimum 
can be treated similarly). Then, we obtain from the strong maximum principle that the 
positive maximum of U w is attained on the boundary T(/), say at zq = (xq, f(xo)). Because 
U m = U w on r(/) the maximum of U m is also attained at zq. Hopf's maximum principle 
ensures then that 



d u U w (zo) < and d u ll m (z ) > 0. 
Moreover, from the last equations in (4.10) we see that 



d v u m 



1 12] 



1 

/' w 



and 



d v u m 



1 



(4.12) 



(4.13) 



Letting x := d u u^ n {zo) and y := d u u 2 n (zQ), we get from the second inequality in (4.12) 
that x > y. Furthermore, since U m attains its maximum at zq its tangential derivative 
d T U m vanishes at zq, meaning that d T u\ n {zo) = d T u 2 n {zo) =: r € R. Whence, |Vu^| 2 (zo) = 



x 2 + r 2 , 



\7u 2 n \ 2 (z ) = y +r 2 and relations (4.13) and (4.12) yield 



< 



fi m {x 2 +r 2 ) fJ, m (y 2 +r 2 )' 



(4.14) 
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But, this is contradicting the fact that the function ]R 3 x i— > x/fi(x 2 + r 2 ) is strictly 
increasing, as we have 

d ( x \ /i m (x 2 + r 2 ) — 2x 2 ji' m {x 2 + r 2 ) 



dx \ Pm (x 2 + r 2 ) J Hm( x2 + r2 ) 

> Vm(x 2 + r 2 ) - 2{x 2 + r 2 )fi' m (x 2 + r 2 ) fr5) 
/^(x 2 +r 2 ) 

where the first inequality holds true only if fj,' m > 0. Thus, U w = and the last equation 
of (4.10) leads us to the desired conclusion F\= F%. □ 

We finish this section with the proof of our main result, Theorem 2.1. The proof follows 
by making use of the abstract theory for fully nonlinear parabolic as presented in Chapter 
8 of [16]. 

Proof of Theorem 2.1. Let a G (0, 1), /i£K, and h be given such that (2.15) is satisfied. 
Then, we infer from Proposition 4.2 that solving the original problem (2.14) is equivalent, 
for small initial and boundary data, to solving the fully nonlinear problem 

dtf = *(*,/), /(0) = /o, (4.15) 

whereby $(t,/) := <Z>(h(t),f) for all (i, /) G [0,T ] x 0, with T > chosen such that 
/i(t) £ W for all t G [0, To]. Since /i is continuous, we obtain from Proposition 4.2 that $ 
and &f are continuous mappings, that is G C([0,To] x 0,Y). Taking into account 

that the above considerations are true for any j3 G (0, a) and that the small Holder spaces 
satisfy the following property 

(/i <J0 (S),^ <J1 (S)) e = HM«»*(S) 

if G (0, 1) and (1 — 0)o~q + 0o~i G M + \N, where (■, -)g denotes the continuous interpolation 
method of DaPrato and Grisvard [5] (see also [16]), the proof of the theorem reduces to 
showing that df&(h, 0) is the generator of a strongly continuous analytic semigroup in 
£(Y). 

In order to prove this property, we infer from Proposition 4.2 that 

d f $(h, 0) = -[d F T(h, 0, 0)]^ 1 o d f T(h, 0, 0), (4.16) 

where dpJ-(h, 0, 0) is the isomorphism given by (4.8). We determine now a representation 
for dfJ-{h, 0, 0). Differentiating (4.2) with respect to /, yields 

d f ?(h,0,0)[f] = d Vm B m (0,h)[d P n0,h)[-7d K (0)[f] + g( Pm - Pw )f]\ for all / G X. 

(4.17) 

Indeed, this simplified expression is a consequence of the fact that T(/)[0, h] = h, TZ(f, h) = 
h, and B m (f,h) = for all / G V. Taking into account that dn(0)[f] = f" for all / G X 
and expending / = ^fcez\{o} o,ke tkx , the relations (4.6) and (4.8) lead us to the following 
representation 

d f $(h t o)f= y, x ^ ikx for /= E a ^ kx ( 4 - 18 ) 

fcez\{o} fcez\{o} 
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A « = : ^ , : :__,2 f1 _M Pi» ~ Pm) ~ k 2 j], k G Z \ {0}. (4.19) 



where the symbol (A&) is given by 

k tanh(fc) 
fj, m (0) + fJ. m tanh 2 (/c) 

The desired generator property may be obtained now from Theorem 3.6.3 in [17] (or The- 
orem 3.4 in [6]). We have verified therewith all the assumptions of Theorem 8.4.1 in [16] 
and the desired statement Theorem 2.1 follows from this abstract result. □ 
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